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Conventional superexchange rules predict ferromagnetic exchange interaction between Ni{II) and 
M (M = Mo^ , , Nb'^). Recent experiments show that in some systems this superexchange 
is antiferromagnetic. To understand this feature, in this paper we develop a microscopic model for 
Ni{II) — M systems and solve it exactly using a valence bond approach. We identify the direct 
exchange coupling, the splitting of the magnetic orbitals and the inter-orbital electron repulsions, 
on the M site as the parameters which control the ground state spin of various clusters of the 
Ni{II) — M system. We present quantum phase diagrams which delineate the high-spin and low- 
spin ground states in the parameter space. We fit the spin gap to a spin Hamiltonian and extract 
the effective exchange constant within the experimentally observed range, for reasonable parameter 
values. We also find a region in the parameter space where an intermediate spin state is the ground 
state. These results indicate that the spin spectrum of the microscopic model cannot be reproduced 
by a simple Heisenberg exchange Hamiltonian. 
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I. INTRODUCTION 

In the expanding field of molecular magnetism, tran- 
sition metal complexes have been the main focus of 
study. Room temperature bulk magnetization has been 
achieved with a hybrid organic radical/ V(II) framework 
Q and also with the Cr(III)/V(II) Prussian blue type 
compound 0,0]. Besides, interesting conjectures such as 
the Haldane conjecture |j| have been verified in molecu- 
lar systems 0. The interest in molecular magnetism has 
expanded in recent times to include low-dimensional ar- 
chitectures exhibiting slow relaxation of magnetization 
known as Single Molecule Magnets (SMMs) [Q], conduct- 
ing magnets and light-triggered magnets The 
large majority of these molecular compounds are based 
on paramagnetic 3d metal ions, the second and third 
row transition metal ions have been envisaged only re- 
cently for the construction of magnetic supramolecular 
compounds . These ions appear however very promis- 
ing in molecular magnetism yielding magnets with fairly 
high ordering temperatures 0|, new SMM's ^3 ^^'^ effi- 
cient systems with light triggered magnetization changes 
[l^. The 4d and 5d ions are characterized by spatially 
more extended valence orbitals, the extension following 
the trend 3c? < 4c? < 5c?. A consequence of this 
ordering of the d-orbital is that the on-site electron re- 
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pulsion is decreased as we go down the column in the 
periodic table. Besides, the metal-ligand bonds become 
more covalent resulting in more efficient electron delocal- 
ization. A series of experimental results obtained on sim- 
ple cyano-bridged heterometallic species formed by self- 
assembling of an octacyanometallate {{Mo^ {CN)s}^~ 
or {W^{CN)s}^^) with a {Nf^L} complex (L = macro- 
cyclic ligand) showed that the nature of the effective su- 
perexchange through the cyano-bridge depends upon the 
actual spin topology of the complex that is being studied. 
For instance, whereas significant ferromagnetic interac- 
tions are found for linear {Ni — M ~ Ni} compounds (M 
= Mo^ or in the [{M{CN)s}] unit), an antiferro- 
magnetic behavior is observed for a cyclic tetranuclear 
compound (Figure^ ^ related larger spin clusters 
of formula [{NiL}i2{Nb{CN^i^e] exhibits an even more 
complex magnetic behavior |l4| . 

In this paper, we explain this puzzling nature of the 
superexchange interaction between two transition metal 
ions A and B, each containing two nearly degenerate or- 
bitals, by developing a microscopic model which natu- 
rally admits antiferromagnetic exchange interactions be- 
tween them. The ion A has two electrons while B has one 
electron; these correspond to the case of Ni^^ for the A 
ion and either Mo^ , or Nb^^ for the B ion. The pa- 
rameters entering the definition of the microscopic model 
are the on-site and inter-site parameters. The on-site pa- 
rameters are both one and two-electron parameters. The 
one-electron parameters are energy differences between 
the active magnetic orbitals on A and B, and the splitting 
of the two orbitals on B 15]. The two electron parame- 
ters are the electrostatic on-site Hubbard repulsion terms 



2 




FIG. 1: Views of the supramolecular compounds top: 
[{NiL}2{M{CN)8}]-^ and bottom : [{NiL'}2{M{CN)s}2]. 



and the inter-orbital Coulomb and exchange integrals. 
The inter-site one-electron parameters are the transfer 
integrals and the inter-site two-electron Coulomb interac- 
tions in the zero differential overlap approximation. We 
identify the critical interaction parameters that control 
the sign of the superexchange interaction and present a 
quantum phase diagram that shows a crossover in the 
ground state from effective ferromagnetic interaction to 
effective antiferromagnetic interaction between the ions. 
We show this behavior in different organizations of the 
AB system such as A — B, linear A — B — A, B ~ A ~ B, 
A - B - A - B, and cycHc A- B - A - B clusters. In 
the next section we introduce the microscopic model for 
superexchange in the Ni^g^^ — M^^^^^ system. This is 
followed by a section in which we describe the valence 
bond (VB) method for solving the model. The last sec- 
tion deals with the results and their discussion. 



II. MODEL HAMILTONIAN FOR 
SUPEREXCHANGE IN AB SYSTEMS 

The nature of exchange interaction between two ions is 
the outcome of a competition between delocalization of 
the electrons which reduces the kinetic energy and direct 
exchange interactions which reduces the electron-electron 



repulsion energy for parallel spin alignment. Usually, 
the kinetic energy is decreased if the electrons are in 
a low-spin state; in the high spin state, the delocaliza- 
tion of electrons is blocked because of the Pauli exclu- 
sion. Besides, low-spin configuration also usually affords 
more phase space for electron delocalization. This is why 
the kinetic exchange is usually antiferromagnetic in na- 
ture. However, when there is a near degeneracy of the 
partially occupied orbitals on an ion, then delocalization 
is possible even when the spins are aligned parallel and 
indeed, the direct exchange interaction favours parallel 
electron spin alignment on such degenerate site. In this 
case, the final outcome cannot be easily predicted and 
depends upon the actual values of the interaction pa- 
rameters. Thus, a microscopic model for explaining ob- 
served sign of the exchange interaction should include 
direct exchange interactions between nearly degenerate 
orbitals on a given site and the relative splitting of the 
degenerate orbitals. Besides these terms, another inter- 
action term of considerable import is the strength of the 
intra-orbital electron repulsions relative to the strength 
of on-site inter-orbital electron repulsions. This is be- 
cause, if the inter-orbital intra-site electron repulsions is 
weak compared to the intra-orbital electron repulsion, it 
will favour single occupancy of two orbitals on the same 
site over double occupancy of one of the orbitals. Oc- 
cupancy of the orbitals controls the nature of the su- 
perexchange process, thus the relative strengths of on- 
site inter and intra orbital repulsions become very impor- 
tant. The other interaction terms which are comparable 
in strength to these repulsion integrals are the intersite 
electron-electron interactions within the zero differential 
overlap approximation and the intra-site electron repul- 
sion between a charge density in one orbital and a charge 
density in the overlap cloud of two orbitals, which we call 
the W term. The model Hamiltonian for investigating 
the superexchange interaction can be written as: 

H = ^ tij{Eij + Eji) 

i <ij> 

—ni{n, - 1) + 

i 

'^^{Ui.vn^hi, -I- -^^^ \{Ei^i, -I- Ei>^i){h., + hi') 
+ + h^,){E,,,, + E,,,,) - 2{[E,^,, + E,, ^i)] 
~^~2~ (-^*^«'-^».''' + Eii^iEii,i + Eis'Ei>,i 
+E.i>^iE^^i, - fii - hi>)} + ^ ^h^hj-, 

= H = AaO'j,^ (1) 

(7 (J 

where, the operators a\ ^ {(ii,a) creates (annihilates) an 
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FIG. 2; Schematic diagram of the active orbitals and the interactions in the A-B-A system. The transfer terms included in the 
model are shown by filled double-headed arrows between pairs of orbitals. Since, only two orbitals are considered on each site, 
the interaction terms are labeled by the site indices (A/B). 



electron in orbital i with spin a. The first line corre- 
sponds to the non-interacting part of the Hamiltonian 
with £i being the energy of the i*'* orbital, and t the 
transfer integral between an orbital on one site and an- 
other on the neighbouring site. All values of t are as- 
sumed to be the same and the orbital energy of the A 
type atoms are fixed at zero, while the orbital energies 
of the B type atoms are —A and (—A -|- 5b) as shown 
in Fig. 121 The second line corresponds to the intra- 
orbital interaction term, with Ui being the Hubbard pa- 
rameter. The remaining lines except the last represent 
the inter-orbital on-site electron repulsion terms, Ui i' = 
W^] = //</':(l)0,(l)^C,(2)0,,(2)d3rid3r2, Ww = 
= / /0*(l)0,(l)4</>*(2)^,,(2)d3rid3r2 and Jw = 

W] = 1 1 r^il)^A^)^, 0*(2)</>,,(2)d3rid3^2, where 
Jul is the exchange integral. Since we deal with only two 
types of ions A and B, and only two orbitals on each 
ion, we label these as Ua, Ua-a, Wa and Ja for the 
A ions and similarly for the B ions. The last line cor- 
responds to the inter-site inter-orbital repulsion integral 
V^, = Nbj] = I I rd^)(^^a)^rJmM2)d^r^^r, and 
is parametrized u sing Ohno parametrization 16] and ap- 
propriate scaling [iTj . 



III. SOLUTION OF THE SUPEREXCHANGE 
MODEL 

The Hamiltonian in eqn. ^is nonrelativistic and hence 
conserves total spin. Since we are interested in the solu- 
tion of the model in different total spin sectors, we use 
the valence bond technique for solving the model. In this 
technique, the complete and linearly independent basis 
states with a given total spin can be generated using ex- 
plicit spin pairings 0. For example, if singly occupied 
sites i and j are spin paired, then a line is drawn between 
sites i and j in the VB diagram to indicate the spin pair- 
ing, {ail3j—Piaj)/y/^, i < j. We say that the line begins 
at site i and ends at site j. If the spins at sites ii,i2 - ■ - k 
are not paired, we pass an arrow through these sites in 



the VB diagram and this denotes the state ai^ ■ ■ - a^. 
Because the Hamiltonian conserves S^^f^i besides S^, it 
is sufficient to work in the subspace Ms — S. In the 
S = subspace, a VB diagram has either empty and 
doubly occupied sites (represented by dots and crosses 
respectively) or lines between singlet paired singly oc- 
cupied sites. A VB diagram involving N orbitals can 
be drawn by arranging the orbitals at the vertices of a 
regular iV-gon and drawing straight lines between ver- 
tices, the electrons at whose orbitals are singlet paired. 
The Rumer-Pauling rule states that all such VB diagrams 
with no intersecting lines form a complete and linearly 
independent set [l9|. The complete and linearly inde- 
pendent set of VB states in nonzero spin space can be 
obtained by taking recourse to modified Rumer-Pauling 
rules [l^ . Some typical VB diagrams are shown in Fig. 
O Each orbital in the VB picture has one of four possi- 
bilities; (i) the orbital is empty, (ii) the orbital is doubly 
occupied, (iii) a line begins at the orbital or (iv) a line 
ends at the orbital. It is possible to associate these four 
possibilities of an orbital in a VB diagram with the four 
states of two bits in an integer, identified with the or- 
bital. Thus, for a eight orbital system, the VB diagrams 
are represented by a sixteen bit integer and these integers 
can be generated in an ascending order. In any given spin 
space, the effect of the operator term Eij in the Hamilto- 
nian on a basis state is to alter the orbital occupancy of 
orbitals i and j (subject to Pauli principle) and pair the 
spins in the orbitals that were involved to yield a new VB 
digram with a known amplitude. If the new VB diagram 
violates Rumer-Pauling rules, it is trivially possible to 
express it as a linear combination of the basis VB states. 

Using the above procedure, the Hamiltonian matrix 
can be set-up in the chosen total spin sector. The re- 
sulting Hamiltonian matrix is sparse (since the number 
of terms in the Hamiltonian is far smaller than the di- 
mension of the complete VB space) and nonsymmetric 
and can be partially diagonalized to obtain a few of the 
low-lying eigenstates, using Rettrups modification of the 
Davidson algorithm, when the Hilbert space dimensional- 
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FIG. 3: Representative VB diagrams in (a) S=0, (b) S=l, 
(c) S=2 and (d) S—3 subspaces. Dots (•) represent empty 
sites and crosses (x) represent doubly occupied sites. In each 
VB diagram, numbers 1 through 8 are the site indices and the 
values at the bottom against each site gives the bit state of 
the bits representing the site in the integer I. The bit-pattern 
and the value of the integer I, representing the VB diagram 
are also shown. 



ity is large |2fl|,|21|. In our case, since the dimensionalities 
of the different subspaces are fairly small (the largest sub- 
space encountered is of dimensionality 1512 for the S—1 
subspace in the A-B-A-B system), full diagonalization of 
the Hamiltonian matrix is resorted to in each spin sector 
to obtain the complete eigenvalue spectrum. 



IV. RESULTS AND DISCUSSION 

The range of values on-site electron repulsion inte- 
grals take is known for the transition metal ions and 
it is reasonable to assume that Ua, when A represents 
Ni^+ is smaller than Ub when B corresponds to M 
{Mo^, W^, Nb^^), since the higher charge on the M 
ion should make the orbitals more compact, resulting in 
larger intra-orbital electron repulsion integrals. It is rea- 
sonable to assume Ua to be 6 eV and Ub to be 8 eV. 
Indeed, we have also verified that the results we present 
do not change qualitatively when these parameters are al- 



lowed to vary by up to 2 eV about the mean values. The 
parameters Ua-a and Ub-b are slightly smaller than the 
corresponding intra-orbital repulsion integral. The inte- 
grals Wa and Wb are much smaller than the Ua~a or 
the Ub-b integrals as they involve overlap charge den- 
sities and both Wa and Wb have been fixed at 1 eV. 
The exchange integrals J a and Jb are somewhat smaller 
than the integrals Wa and Wb- We also note that since 
it is experimentally known that the spin on the A ion 
{Ni^^) is always one we fix J a at a somewhat large value 
of 0.7 eV. Large J a reduces the repulsion between elec- 
trons with parallel spin alignment compared to antipar- 
allel spin alignment in the two orbitals on site A leading 
to a spin-1 object on site A and this holds true even 
when the degeneracy of the two orbitals on the A site is 
slightly lifted. For this reason, we have assumed the two 
orbitals on A to be degenerate. In each case, we have 
also verified that the expectation value of S\ operator is 
nearly 2.0 confirming that the spin on the A site is very 
nearly one in all cases. The total spin of the system in 
the ground state is sensitive to the parameters 5b, Jb 
and Ub-b- A large 5b would result in an antiferromag- 
netic exchange interaction since the virtual state with a 
doubly occupied lower orbital on the B site has a lower 
energy than spin one state on the B ion. A large Jb 
would however lower the energy of the virtual state in 
which the B ion has a spin one configuration. A small 
Ub-b would also favour a spin one virtual state on the 
B ion by favouring single occupancy of the two active B 
ion orbitals. We have solved the model Hamiltonian over 
a wide range of the parameters and have obtained the 
quantum phase diagrams for demarcating the low-spin 
and high-spin ground states in this parameter space. 

In Figures 01 |H and El we have presented the contours 
of superexchange values (Js) in these systems for vari- 
ous values of the model parameters obtained by fitting 
the low-spin - high-spin gap to a spin Heisenberg Hamil- 
tonian involving spin-l's at A sites and spin-l/2's at B 
sites. The solid line corresponding to Jg = provides 
the phase boundary between the high-spin and low-spin 
ground states. The contours corresponding to fixed Js 
values and represented by dotted lines are obtained by 
spline interpolation using MATLAB. We note that small 
Ub-Bj small 5b and large Jb values promote a high spin 
ground state while large Ub-b, large 5b and small Jb 
values promote a low spin ground state. We also note 
that the superexchange Js values, for the same model 
parameters, are larger for smaller system sizes. 



A. Phase diagram in the Ub-b — Jb plane 

In the Ub-b — Jb plane, we should expect to see 
the high-spin ground state for large Jb and small Ub-b- 
This is because the large exchange integral on the B site 
will favour parallel alignment of the electrons on the B 
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FIG. 4: Contours of the effective superexchange constants (Js cm^^) of (a) A — B, (b) A — B — A, (c) B — A — B chains and 
(d) cyclic A — B — A ~ B systems as a function of Ub-b and Jb- The phase diagrams are obtained for t=0.1 eV, A=0.0 eV; 
i5b=0.0 eV; Ua^Q eV; Ub=9> eV; Ua-a=^ eV; Ja=0.7 eV; Wa^Wb=1 eV. All the systems display high-spin ground state at 
higher Jb and lower Ub-b values. 



site, if the virtual transfer of the electron does not lead to 
a doubly occupied B site orbital. The latter requirement 
is guaranteed if Ub-b is small. In Fig. ^we see that in 
all the cases, the ground state corresponds to the high 
spin state for small Ub-b and large Jb- We also note 
from the shift in the phase line to the right, that with 
increase in system size the high-spin state is favoured 
for larger Ub-b values. From the gap between the high 
spin and the low spin states, we have also calculated 
the effective superexchange parameter. We find that the 
largest ferromagnetic superexchange J value is 101.42 K 
(70.49 cm-i) for Ub^b = 5 eV and Jb = 0.5 eV, while 
the largest antiferromagnetic superexchange J value cor- 
responds to -33.36 K (-23.19 cm"!) for Ub-b = 7 eV 
and Jb =0.1 eV in the linear A — B cluster. Another 
interesting feature to note is that in the B — A~ B clus- 
ter, the low-spin ground state appears for a much smaller 
Ub-b value, for a given Jb than in the A — B — A cluster. 



B. Phase diagram in the Ub-b — 5b plane 

The phase diagram in Ub-b — <5_b plane for all the 
four systems is shown in Fig. [S] We note that the ground 
state spin is extremely sensitive to the 5b value. Even 
a small splitting of the B site orbitals forces the system 
into a low-spin ground state. This is because the lifting 
of the degeneracy of the B site orbital favours doubly 
occupied lower energy orbital in the virtual state, which 
would result in stabilization of the low-spin state. Here 
again, we note that the high spin ground state is more 
robust in the larger clusters. In the cyclic A— B — A — B 
cluster, for Ub-b — 5.5 eV, the ground state shifts to 
low-spin state for Sb = 0.063 eV while in the A — B — 
A cluster, this occurs at Sb = 0.04 eV and in the A — 
B system the ground state ceases to be the high spin 
state for 5b > 0.016 eV. The highest ferromagnetic 
superexchange constant J is observed for 5b = and 
Ub-b = 5 eV while the highest antiferromagnetic J is 
observed for 5b = 0.1 eV and Ub-b = 7 eV in all the 
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FIG. 5: Contours of the effective superexchange constants {Js cm^^) of (a) A — B, (b) A — B — A, (c) B ~ A ~ B chains and 
(d) cyclic A ~ B — A — B systems as a function of Sb and Ub-b- The phase diagrams are obtained for t=Q.l eV, A=0.0 eV; 
Ua=6 eV; Ub=8 eV; C/a-a=4 eV; Ja~0-7 eV; Js=0.3 eV; Wa=Wb=1 eV. All the systems display high-spin ground state at 
lower 5b and Ub-b values. 



clusters. We also can note that for a given 5b value, we 
have a low spin ground state in the B — A — B system, at 
much smaller value of Ub-b value than in the A — B — A 
system. 



C. Phase diagram in the Jb Sb plane 

Large values oi Sb has the effect of promoting the an- 
tiferromagnetic superexchange while a large Jb favours 
ferromagnetic superexchange. Thus, we see from Fig. |H| 
that the high spin state is the ground state below the 
phase line while the low-spin state is the ground state 
above the phase line. The phase line shifts higher in the 
Jb — Sb plane as the system size is increased thereby 
showing that the high spin ground state is more stable 
to lifting of the degeneracy of the orbitals on the B site 
for larger system sizes. In this plane, the behaviour of 
the A — B — A and the B — A — B systems are almost 
identical. 



One very interesting feature not clearly seen in the 
phase diagrams of the cyclic A — B — A — B cluster is the 
appearance of the intermediate spin S=2 ground state 
over a very narrow range of parameter values between 
the S=l and the S=3 ground states. The intermediate 
spin ground state is seen only in the largest cluster we 
have studied. In fact, in the linear A~ B ~ A— B cluster, 
this region extends over a wider parameter values, as 
seen in Fig [T] Thus, it may be possible to synthesize 
high nuclearity complexes in intermediate spin ground 
states. The physics behind the existence of such ground 
states is however somewhat different from the frustrated 
magnetic exchanges present in SMMs. In our case, it is 
difficult to identify the exchange interaction between two 
sites as either ferro or antiferromagnetic and magnetism 
can be viewed as a whole and not pair- wise, as is usually 
the case. Besides, it is also not possible to map our model 
onto a simple Heisenberg exchange Hamiltonian. 

To conclude, we have developed a model Hamiltonian 
which admits both low-spin and high-spin ground states 
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FIG. 6: Contours of the effective superexdiange constants (Js cm~^) of (a) A — B, (b) A — B — A, (c) B — A — B cliains and 
(d) cyclic A — B — A — B systems as a function of Jb and 5b- Tlie phase diagrams are obtained for t=0.1 eV, A=0.0 eV; 
Ua=6 eV; f/s=8 eV; (7a-a=4 eV; (/b_s=5.5 eV; Ja=0.7 eV; Wa=Wb = 1 eV. AU the systems display high-spin ground state 
at lower 5b and higher Jb values. 



for small changes in the values of the model parame- 
ters. The only model parameter that can perhaps be 
determined directly for a system is (5s, from spectro- 
scopic data. Other parameters can only be inferred indi- 
rectly either from electron spectroscopic studies on sim- 
pler systems or from abinitio model system calculations. 
Our model can explain the observed antiferromagnetic 
exchange in the NbfiNii? and related systems, contrary 
to superexchange rules 0, 0| . The model also yields 
reasonable effective superexchange constants for model 
parameters in the accepted range. 
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